Dynamics of Current Induced Magnetic Superstructures in Exchange- Spring Devices 
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Thermoelectric manipulation of the magnetization of a magnetic layered stack in which a low- 
Curie temperature magnet is sandwiched between two strong magnets (exchange spring device) is 
considered. Controllable Joule heating produced by a current flowing in the plane of the magnetic 
stack (CIP configuration) induces a spatial magnetic and thermal structure along the current flow 
- a magneto-thermal-electric domain (soliton). We show that such a structure can experience 
oscillatory in time dynamics if the magnetic stack is incorporated into an electric circuit in series 
with an inductor. The excitation of these magneto-thermionic oscillations follow the scenario either 
of "soft" of "hard" instability: in the latter case oscillations arise if the initial perturbation is large 
enough. The frequency of the temporal oscillations is of the order of 10 5 -7- 10's _1 for current 
densities j ~ 10 6 10 7 A/cm 3 . 

PACS numbers: 



I. INTRODUCTION 

The electric control of magnetization on the nanome- 
ter length scale became a subject of intensive study af- 
ter the seminal suggestion by Slonzewsky 1 ,2 to use the 
spin torque transfer (STT) technique, where a current of 
spin-polarized electrons is injected into a magnetic ma- 
terial. Even though the magnetic tuning induced by the 
exchange interaction between the injected electrons and 
those that make the material magnetically ordered can 
only be achieved in a very small region, large current 
densities are needed to get a significant STT effect. Such 
current densities can be reached in electric point contacts 
where values of the order of 10 8 4- 10 9 A/cm 3 have been 
achieved and the STT effect was detected. 3,4 By a further 
increase of the current thermal Joule heating of the ma- 
terial becomes unavoidable. However, the fact that such 
heating can be precisely controlled by the bias voltage 
allows thermoelectrical manipulation of the magnetiza- 
tion direction based on the orientational phase transition 
in layered magnetic structures predicted in Ref. 5-7. It 
was shown there that in such structures with different 
Curie temperatures of the layers there can exist a finite 
temperature interval inside which the angle of the relative 
orientation of the layer magnetization 9 depends on tem- 
perature, 8 = Q(T), changing reversibly with a change of 
temperature in the whole interval of angles from the par- 
allel to the antiparallel orientation. As a result, by con- 
trolling the Joule heating by the bias voltage V one may 
smoothly and reversibly control the relative magnetiza- 
tion angle in the whole angle interval < 0[T(V)] < ir. 
There was also suggested a thermal-electronic oscillator 
device based on this magneto-thermal effect in magnetic 
exchange spring structures. 

In paper 8, for exchange-spring magnetic stacks'* with 



an in-plane electric current flow (CIP) (see Fig. 1), it was 
predicted and investigated a possibility of generation of 
magnetic spatial superstructures (domains and solitons) 
the appearance of which and their parameters can be con- 
trolled by the bias voltage. In the present paper we study 
temporal dynamics of such spatial structures in the case 
that the exchange-spring magnetic stack is incorporated 
in an external circuit in series with an inductor 10 . Such 
oscillations can be understood as a realization of the bi- 
stability of an N-shaped IV curve of the device in the 
regime at which a large enough inductance of the circuit 
determines slow temporal variations of the current. 

The structure of the paper is as follows. In Section II 
we present in short some peculiar features of the depen- 
dence of the magnetization orientation angle in the stack 
on temperature T and some properties of the magnetic 
stack with a magnetic-thermal-electric domain (MTED) 
inside it needed for further considerations. In Section 
III we show that MTED loses its stability if the stack is 
incorporated in an external circuit with a large enough 
inductance. We develop there an adiabatic theory which 
allows an analytical description of the MTED dynamics. 



II. THE MAGNETIC-THERMAL-ELECTRIC 
DOMAIN IN THE STACK UNDER JOULE 
HEATING 

The system under consideration has three ferromag- 
netic layers in which two strongly ferromagnetic layers 
and 2 are exchange coupled through a weakly ferromag- 
netic spacer (layer 1) as is shown in Fig. 1. We assume 

that the Curie temperature Tc of layer 1 is lower than 
the Curie temperatures Ti°' 2 ' of layers 0, 2; we also as- 
sume the magnetization direction of layer to be fixed; 
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FIG. 1: Orientation of the magnetic moments in a stack of 
three ferromagnetic layers (0, 1, 2); the magnetic moments 
in layers 0, 1, and 2 (shown with short arrows) are coupled 
by the exchange interaction thus building an exchange spring 
tri-layer; H (shown with the arrow outside the stack) is an 
external magnetic field directed antiparallel to the magneti- 
zation in layer 0. The current J (shown with the long arrow) 
flows in the layer plane (that is along the x-axis). 



this stack is under an external magnetic field H directed 
opposite to the magnetization of layer 0. We require this 
magneto-static field to be weak enough so that at low 
temperatures T the magnetization of layer 2 is kept par- 
allel to the magnetization of layer due to the exchange 
interaction between them via layer 1. At H = and 
T > T^ this tri-layer is similar to the spin-flip "free 
layer" widely used in memory device application J 1 . The 
stack is incorporated into an external circuit along which 
the total current J flows through the cross-section of the 
layers: 
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where R(0) and i?o are the magnetoresistance and the 
angle-independent resistance of the stack, is the angle 
between the magnetization directions of layers and 2, 
V is the voltage drop across the stack. 

In paper 1 ', it was shown that parallel orientations of the 
magnetization in layers 0, 1 and 2 becomes unstable if the 
temperature exceeds some critical temperature Tc < 
Tc , The magnetization direction in layer 2 smoothly 
tilts with an increase of the stack temperature T in the 
temperature interval Tj™' < T < Tc . The dependence 
of the equilibrium tilt angle between the magnetization 
directions of layers and 2 on T and the magnetic filed 
H is determined by the equation'' 



6 = D(H 7 T) sinG, T<T C (1) 

e =±tt, t>t c (1) 
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Here L Xl Mi(T) and L 2 , M 2 (T) are the widths and 
the magnetic moments of layers 1 and 2, respectively; 
ol\ ~ /i/aMf(0) is the exchange constant, 1\ is the ex- 
change energy in layer 1, [i^ is the Bohr magneton, Ub 
is Boltzmann's constant, and a is the lattice spacing. 
D(H,T) is a dimensionless parameter that determines 
the efficiency of the external magnetic in the misorienta- 
tion effect under the consideration. It is a ratio between 
the energy of magnetic layer 2 in the external magnetic 
field and the energy of the indirect exchange between 
layers and 2 (see Fig. 1). 

The critical temperature of the orientation transition 12 
Tc° r ^ is determined by the condition D (T) = 1 and is 
equal to 
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Taking experimental values (see Ref. 6) Li = 30 nm, 
L 2 = 12 nm, T c (1) = 373K, and magnetic field H = 10 -f- 
47 Oe one finds Dq = 0.1 0.36 and the temperature 
interval ST = T c (1) - T c (or) = D T^ ] w 37.3 ~ 134 K. 

If the stack is Joule heated by current J its tempera- 
ture T(V) is determined by the heat-balance condition 



JV = Q(T), J = V/R eS {Q), 



where 



R{Q)R a 
R(Q) + R ' 



(6) 



(7) 



and Eq. (2) which determines the temperature depen- 
dence of 0[T]. Here V is the voltage drop across the 
stack, Q(T) is the heat flux flowing from the stack and 
Refi(®) is the total stack magnetoresistance. Here and 
below we neglect the explicit dependence of the magne- 
toresistance on T considering the main mechanism of the 
stack resistance to be the elastic scattering of electrons 
by impurities. Below we drop the subscript "eff" at the 
magnetoresistance symbol R e g. On the other hand, we 
consider the temperature changes caused by the Joule 
heating only in a narrow vicinity of Tc which is suffi- 
ciently lower than both the critical temperatures Tj ' 2 ^ 
and the Debye temperature. 

Equations (6) and (2) define the IVC of the stack 
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where Q(V) = Q[T(V)]. The differential conductance 
Gdiff = dJ/dV is negative 6 if 
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FIG. 2: Current-voltage characteristics of the magnetic stack 
of Fig. 1 in which the magnetization directions in the layers 
are homogeneously distributed along the x— direction (that is 
along the stack). It was calculated for R(0) — R+—R- cos 0, 
R-/R+ = 0.2, Do = 0.2; J c = V c /R(w). The branches - a 
and b — b' of the IVC correspond to parallel and antiparal- 
lel orientations of the stack magnetization, respectively (the 
parts a — a' and — b are unstable); the branch a — b corre- 
sponds to < Q[T(V)] < ty. 



where D = (T/Q)(dQ/dT)D « D In this case the 
current- voltage characteristics of the stack (IVC) is re- 
shaped as is shown in Fig. 2. 

In paper Ref. 8 it was shown that a homogeneous in 
space distribution of the magnetization direction, tem- 
perature and electric field along the spring-type magnetic 
stack becomes unstable and a magneto-thermal-electric 
domain (MTED) spontaneous arise 14 if the electric cur- 
rent flows in the plane of the layers (CIP-cofiguration) 
and the length of the stack exceeds the critical length 




(10) 



where K is the thermal conductivity, p(0) is the magne- 
toresistivity and r = (p(0) — p(iv))/(p(0) + p(ir)). Pa- 
rameters of the MTED which arises in the stack, de- 
pend on the stack length and the voltage drop across it 
that allows to create and control magnetic structures in 
magnetic devices. Using Eq. (10) and the Lorentz ration 
k/ct = Tr 2 kgT/3e (a = p^ 1 and e is the electron charge) 
one finds that 
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for a realistic experimental situation- 1 ' r ~ 0.1 -r 0.3, 
TP ~ 10 2 K, p(0) ~ 10~ 5 fh:m, j ~ 10 6 10 7 A/cm 2 . 

The magnetic-thermal-electric domain (MTED) which 
determines the space distributions of the temperature 



Td(x) and the magnetization direction angle <d(x) which 
spontaneously arise in the magnetic stack, satisfies the 
following equation s : 



dx 2 



j 2 p{6(T d )} - Q(T d )/n st 



in which current j is found from the equation 
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j (p[Q(T d )}) = j 
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supplemented with the periodic boundary condition 
Td{x + L) = Td(x) (x is the coordinate along the stack, 
L and Sl s t are the length and the volume of the stack, re- 
spectively; the brackets < ... > indicate an average over 
x along the whole stack of the length L). 

The space distribution of the temperature Td(x) and 
the magnetization direction angle Q(x) in a long stack 
(L ^> L c ) with a MTED inside it is shown in Fig. 3. 

In this limit (L > L c ) the MTED is of a trapezoid 
form with two planar segments of the length Lj and L/j 
(Li + Ln = L) with two transition regions of the width 
- L c (see Fig. 3). 

Neglecting the contribution of the transition region in 
the total voltage drop across the stack, the dependence 
of the lengths Lj and Ljj on the voltage drop V across 
the stack can be written as 



V/L-Si _ Sn-V/L 

L I ~ — n ~ J L H ~ 
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where £/ = p(0)j, £n = p(n)j and r = p~/p+, p± = 
(p{tt) ± p(0))/2. In this approximation these formulas 
are valid in the range of V where Ljjj > 0. The current 
J = Sj and the voltage V = £L arc coupled via current- 
voltage characteristics j = jd{V/L) for the stack with 
MTED inside it. This dynamic IVC has a form close 
to a plateau in which the current j coincides with the 
stabilization current s jo to an accuracy exponential in 
L/L c . The dynamic CVC can be presented in an implicit 
form as 

£ - P(0)j L 



j - jo = cxp ^ — ^-r\ 

(£- p(n)j L \ 

-Juexpi— r — \ 

1 r {JoP+) L a ) 



(15) 



Here £ = V/L, Lq = v/ f' T {Tc )/n ~ L c and the con- 
stants Jiji are of the order of j while fx = df/dT. 



III. TIME EVOLUTION OF 
MAGNETO-ELECTRIC-THERMAL DOMAIN IN 
MAGNETIC STACK 

Here we assume that the magnetic stack is placed in 
series with an inductor. We also assume the bias- voltage 
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FIG. 3: The coordinate dependence of the temperature, 
T(x), and the magnetization direction, Q(x), in the mag- 
netic stack with a magneto-thermal-electric domain inside it, 
Lu is the length of the "hot" part of the MTED which is 
defined in such a way that the length of the "cold" part is 
Li — L — Lu Calculated are made for R(Q) — R+ — R- cos©, 
R-/R+ = 0.2, Do = 0.2; J/J c = 1.0265, J c = V c /R(tv), 

L = 



y/j*p(n)/KT™ 



regime that is the resistance of the external circuit in 
which the magnetic stack is incorporated can be ne- 
glected in comparison with that of the stack. In this 
case, taking into account that the temperature T{x,t) 
(here t is the time) satisfies the equation of continuity of 
heat flow (see, e.g., Ref. 29), one obtains the following 
set of basic equations of the problem: 



dT d / 0T\ 
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+ (p[Q(T)])j(t) = j 



f(T,j) = Q(T)/n st -j 2 (t)p{0(T)] 



(16) 



(17) 



Here the dependence &(T) is given by Eq. (2). In the 
above equation j(t) = J(t)/S is the current density which 
is independent of x due to the condition of local electrical 
neutrality, 5 is the cross-section area of the stack and L is 
its length, c v is the heat capacity per unit volume, Q(T) 
is the heat flux flowing from the stack, Q, st is its volume 
and p[8] = R[9]S/L is the stack magneto- resistivity (see 
Eq. (7)), C is the inductance, and V is the bias voltage; 
the definition of the brackets <...> is the same as in 
Eq.13. 

The boundary condition to Eq. (16) is the continuity 
of the heat flux at the both ends of the stack (which is 
coupled to an external circuit with a fixed voltage drop V 
over the stack). We shall not write down the expression 
for this condition, since the domain dynamics does not 
sufficiently depend on it. Instead, for the sake of simplic- 
ity we use the periodicity condition T(x + L,t) = T(x, t). 



The set of equations Eq. (16) always has a steady-state 
solution which represent a magnctic-thcrmal-clectric do- 
main Td{x) which satisfies Eq. (12). A study of the sta- 
bility of this solution carried out using a linearized sys- 
tem Eq. (16) shows that the MTED loses its stability 
if the inductance C exceeds some critical value C c \ ~ 
ttR+(L/L c ) expL/L c where tt ~ CyT jRJ 2 is the char- 
acteristic evolution time of the temperature . Therefore, 
for the case L ^> L c this instability develops under condi- 
tions that the characteristic time of the current evolution 
T c ~ C/R is much longer than the temperature evolution 
time tt (their ratio is tc/tt ~ (L/L c )expL/L c . Then 
any small fluctuation develops in such a way that after 
a time t ^> tt the domain recovers its initial trapezoidal 
form and only its length Lu varies slowly with time t. 
This fact makes it possible to obtain a reduced descrip- 
tion of the nonlinear dynamics of a MTED as follows. 
We seek the solution of the first equation in Eq. (16) in 
the form 



T{x,t) =T d {x l3d {t))+T 1 {x,t) 



(18) 



where T^(x,jd) is the solution of the equation which is 
obtained from Eq. (12) by the substitution j — > jd(t); the 
parameter jd governing the length of the domain Lu(jd) 
is a slow function of the time t which has to be found; 
Ti(x, t) is a correction which is small with respect to T d . 
Inserting Eq. (18) into the first equation in Eq. (16) and 
linearizing it one gets the following equation for T\\ 



HTi 



where H is a Hermitian operator with a periodic bound- 
ary condition: 



H = -n— + ^ - j 2 d p T (T d ) + Q T (T d )/n st j (20) 

(here p{T) = p[<d(T)} and F' T = dF/dT). The require- 
ment of the smallness of the ratio \Ti/T d \ <C 1 allows to 
find the needed equations describing the adiabatic MTED 
evolution. 

The solution of Eq. (19) is found in the form of expan- 
sion 



T 1 (x,t) = ^2A u (t)9 v (a 



(21) 



where ^ v are the eigenf unctions of H which satisfy the 
Sturm-Liouville equation 



K ¥x~ ~ 3dpT ^ Td) + QT(Td)/^»t\ *v = K^v (22) 



Multiplying the both sides of Eq. (19) by from the 
left and averaging over the period L one finds 

KA V = (( - c y g f + (j 2 - f d ) P (T d ))^) (23) 



where dTd/djd in terms of eigenfunctions is presented 
in Appendix A, Eq. (A9). 

As the eigenvalue Ao is exponentially small (see 
Eq. (A5)), from Eq. (19) it follows that the requirement 
\Ti/T d \ <C 1 is satisfied only when the right-hand side of 
the equation with v = is equal to zero: 



♦ L/L 
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(24) 



Using Eq. (A3) one sees that the right-hand side of the 
equation with v = 1 is equal to zero and hence the factor 
Ai = in the left-hand side of it docs not violate the 
above-mentioned requirement. 

Taking into account the second equation in Eq. (16) 
and Eq. (24) together with Eqs.(A8, All) one finds a set 
of equations which describes the temporal dynamics of a 
MTED: 



d Jd i ■ ^ -3d 
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= V/L-(p(T d (x,j d ))j 



(25) 



Here £ = S £/L and ujq = — Ao/cy (here Ao is defined 
by Eq. (All)) may be written as 
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(26) 



(27) 



being of the order of the characteristic evolution time of 
the temperature tt, and £ d is defined in Eq.(A6). 

Throughout the range of the existence of a trapezoidal 
MTED one has \j — j 1 <C jo and the set of equations 
Eq. (25) can be reduced to an equation for the electric 
field which is coupled to the current j d via the voltage- 
current characteristic £ — £ d (j d ) (see Eqs.(13,A6)). Dif- 
ferentiating the first equation in Eq. (25) multiplied by 
wo with respect to t, and inserting in the resulting ex- 
pression the above-mentioned electric field £ as a new 
variable together with the second equation in Eq. (25) 
one obtains the following equation in the form of a linear 
oscillator with a nonlinear nonconservative term: 



£ 



£ 

£jo 



to d£ 



)e + £L(£- V/L) = 0. (28) 
Cto 



Here j d (£) is the CVC of the magnetic stack with a 
MTED. 

The static points of Eq. (28) and the second equation in 
Eq. (25) correspond to the domain solution £ = V/L and 
j = ]d{V/L). As dj d /d£ < 0, in the range of inductance 
£ > C~ci the factor in front of £ is negative and a domain 
is absolutely unstable. Here the critical inductance is 



£ c \ ~ -t ' — I 



V d£ I £=£i. 



t oP+ cxp(L/L c ) (29) 
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FIG. 4: Temporal evolution of the length of MTED Lu (i) for 
the case £ c i > £ > £ C 2- The unstable limiting cycle (shown 
with a thick line) separate the initial the phase plane into two 
regions: any initial state inside the limiting cycle develops 
into the length of the steady stationary MTED (shown with 
a dot) while an initial outside it results in oscillations of the 
MTED length with an increasing in time amplitude until the 
MTED disappears, that is Lu reaches either Lu » or Lu w 
L. Calculations are made for R-/R+ — 0.2, Do = 0.2 and 
tt/tc = 0.1 where tt and rc are the characteristic times of 
the temperature and current developments, respectively; Jo 
is the stabilization current. 



where £i n f is the electric field at which —dj d /d£ has a 
maximum that is at which dj^/d£ 2 = 0. In this region 
the voltage drop across the magnetic stack V(t) = £(t)L 
oscillates and hence the length of the domain Lu(t) (see 
Eq. (14)) oscillates with an amplitude that increases in 
time. At the time when Lu reaches either values 0, L (de- 
pending on the initial fluctuation) the domain disappears 
and the sample becomes homogeneous with temperatures 
T min or T max , respectively. Further temporal evolution 
of the system is described by the set of equations Eq. (16) 
an analysis of which for the case tt/tc <C 1 shows that 
the system exhibits stable large-amplitude spontaneous 
oscillations of the temperature T, current current J, mag- 
netization direction 0(T), and the voltage drop across 
the stack V = Ji?[0(T)], the oscillations being the same 
as those predicted in Ref. 6 for a magnetic stack in the 
absence of MTED. 

Investigation of stability of the domain solution of 
Eq. (28) with the Poincare method (see, e.g., Ref. 32) 
shows that at SC = C c i — £ <C £ there is an unsta- 
ble limiting cycle in the phase plane (£,£) the radius of 
which K (x y/ S£/£ c i and hence it increases with a de- 
crease of the inductance. From here and from the fact 
that at £ = the domain is absolutely stable follows an 
existence of a second critical value of the inductance £ c i 
which limits the interval of the hard excitation of oscil- 
lations. Therefore, the range of values of the inductance 
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FIG. 5: Temporal evolution of the current J(t) and volt- 
age drop V(t) in a magnetic stack with a magnetic-thermal- 
electric domain (MTED) inside it for the case C c i > C > C C 2- 
The steady CVC of a homogeneous stack CVC with a MTED 
are shown with a dashed and dashed-dotted lines, respectively. 
The stack is in a bistable state: depending on the initial condi- 
tions the system goes either to MTED a stable steady MTED 
(which is shown with a dot) or goes to a stable limiting cycle 
(the largest closed curve) corresponding to spontaneous os- 
cillations of the current J(t), voltage drop V(t), temperature 
T(t), and the magnetization direction O(t), the stack being in 
a homogeneous state (in which the MTED has disappeared). 
There is an unstable limiting cycle that separates the initial 
states which develop either to the MTED or the the oscilla- 
tions as is shown by two black arrows. Calculations are made 
for R-/R+ = 0.2, Do = 0.2 and t t /t c = 0.1 where r T and 
T£ are the characteristic times of the temperature and current 
developments, respectively; J c = V c /R(ir). 



C c i > C > C C 2 is the region in which a stable MTED, an 
unstable limiting cycle, and a stable limiting cycle coex- 
ist. For this case, temporal developments of the MTED 
length La, current J and voltage drop V across the mag- 
netic stack are shown in Fig. 4 and Fig. 5. 

IV. CONCLUSION 

We have considered the temporal evolution of 
magneto-thermo-electric domains which spontaneously 
arise in an exchange-spring magnetic stack with the elec- 
trical current flowing along the layers (CIP configura- 
tion). For the case that such a stack is incorporated in 
an electric circuit in series with an inductor we have found 
critical values of the inductance at which the magneto- 
thermal domain loses its stability. We have developed a 
perturbation theory in an adiabatic parameter (which is 
the ratio between the temperature and current character- 
istic evolution times) theory which allows to investigate 
the evolution of this instability into an oscillatory in time 
regime. The excitation of such magneto-thermionic os- 



cillations follow the scenario of cither "soft" of a "hard" 
instability; in the former case any perturbation results in 
a spontaneous transition to the oscillatory regime while 
in the latter case oscillations only appear if the initial 
perturbation is large enough. The frequency of the tem- 
poral oscillations of the magnetization direction in the 
stack is of the order of 10 5 10 7 Hz. 
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Appendix A 

According to Ref. 8 Td(x) can be written as 

T d {x) = 0(s + x+) + $(x+ - a?) - T max (Al) 

where the function $(x) is a domain- wall type solution of 
Eq. (12) at j = j Q and L —> oo which has the following 
asymptotic 

lim ti(x) = T min 

lim 0(x) = T max (A2) 

x— >oo 

Here ±x+ are the points of deflections of the curve 
dT c i{x)/x so that Lu = x + — (— x+) = 2x + is ap- 
proximately the length of the "hot" section of a trape- 
zoidal MTED having the maximal temperature T max and 
Li = L — Lu is the length of its "cold" section having 
the minimal temperature T m .;„ (see Fig. 3). 
As one sees from Eq. (22) the function 

U(x) = -f dP ' T (T d ) + Q T (T d )/n st 

in the former equation consists of two symmetrical "po- 
tential" wells of the width ~ L c , separated by a barrier of 
the width of order of Lu. In this case (see, e.g., Ref. 31) 
tunneling between the wells results in a splitting of each 
of the eigenvalues of a completely separated well into 
two neighboring ones, the splitting being proportional to 
cxp —Ln/L c , while the corresponding cigenfunctions are 
symmetric and antisymmetric combinations of the corre- 
sponding cigenfunctions of the separated left and right 
wells. 

From the translation symmetry of the time indepen- 
dent equation Eq. (12) it follows that the eigenfunc- 
tions of the Hermitian operator H include the function 
"Fi = Td/dx corresponding to the eigenvalue Ai = (it 
is easy to check inserting this function in the Eq. (22)). 
As it follows from Eq. (Al) the eigenfunction 'Fi(x) is 

dx dx dx 
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where two functions in the right-hand side are eigen- 
functions of the left and right wells when the tunnel- 
ing is ignored. This is an antisymmetric eigcnfunc- 
tion ('4 r (— x) = —ty(x)) corresponding to the eigenvalue 
Ai = and hence the nearest eigenvalue Ao is negative 
and the corresponding cigcnfunction is symmetric: 

Mx) = *Kx + x ± ) d*{x+ - x) (A4) 



(l.T 



dx 



while 



Ao = -K 



( dd{x + x + ) d 2 i9(x + x+) 



\ dx 
oc exp (— L/L c ) 



dx 2 



x=0 



(A5) 



To express the eigenvalue Ao in terms of the differential 
resistivity we use the following reasoning. 

According to Eq. (13) the voltage-current characteris- 
tic of the stack with a MTED in it is 



£ d (j) = j(p[e{T d (j))}) 



(A6) 



Differentiating the both sides of Eq. (A6) and Eq. (12) 
with respect to j one finds the differential resistivity of 
the magnetic stack with a MTED as 



(A7) 



and an equation for dTd/dj 



(A8) 



where the operator H is given by Eq. (20). 

Expanding dTd/dj in the form of the eigenfunctions 
(see Eq. (22)) one finds 



dj 



A„ 



(A9) 



Inserting it in Eq. (A7) one finds 
d£d 



dj 



+ 3*3 £ (Pt(Wp(W (Aiq) 

At/ 



Remembering that Ao is exponentially small and hence it 
gives the main contribution to the sum with respect to v 
one finds 

_ I6r(j 2 p + ) p[Q(T)]dT 
L 2 < (dT d /dx) 2 > d£ d /dj ( ' 
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